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REFRACTION  EQUATIONS  FOR  HYDRONS  CONSIDERING 


The  main  purpose  in  this  paper  is  to  derive  the  refraction  equations  for 
surface  water  waves  which  are  moving  in  a  variable  current  over  a  variable 
bottom  topography.  A  second  purpose  in  the  paper  is  to  compare  the  refrac¬ 
tion  equations  obtained  by  various  investigators.  Although  the  results  are 
presented  for  the  example  of  water  waves,  the  refraction  expressions  obtained 
are  applicable  to  all  kinds  of  waves. 

The  waves  are  presumed  to  propagate  in  the  form  of  wave  packets.  A 
wave  packet  is  obtained  by  the  summation  of  all  sine  waves  over  a  narrow 
range  of  both  frequencies  and  directions.  Synge  (1962)  has  suggested  the 
name  hydron  for  wave  packets  associated  with  water  waves.  It  is  assumed  that 
wave  properties  such  as  the  wave  speed,  direction,  and  wavelength  do  not 
change  much  over  distances  of  one  wavelength.  Ray  curvature  expressions  are 


derived  following  Che  procedure  outlined  by  Laudau  and  Lifshltz  (1959).  The 
water  wave  velocities  are  assumed  to  be  much  greater  in  magnitude  than  the 
currents,  and  terms  higher  than  first  order  in  the  current  velocity  are 
ignored. 

2.  Hydrons  in  currents 


An  unprimed  coordinate  system  xy  is  fixed  in  space  on  the  water  surface. 

A  primed  coordinate  system  x'y'  is  considered  to  be  stationary  with  respect 
to  a  water  current  and  is  moving  with  velocity  u  relative  to  the  xy  coordinate 

system.  The  position  vector  r  in  the  xy  coordinate  system  and  the  position 

*# 

vector  r'  in  the  x'y'  coordinate  system  are  related  by 

7L  m  at  +  n‘  (1) 

where  t  is  time. 


For  an  observer  moving  with  the  current  a  hydron  is  expressed  by 


(Breeding,  1978)  XI  +< 

i  s  ui-a  -<•>■*)  (  ft 


(  -  cr't) 

\  A(w)  e 


n-t 


where  If  is  the  displacement,  A  is  the  amplitude,  k  is  the  wave  number,  and 
(j)  is  the  radian  frequency  of  the  wavelets.  The  wave  number  m  and  the  radian 
frequency  O’  of  the  wave  packet  are  defined  by 

wv.  - 

O'  =  AU> 

The  average  radian  frequency  of  the  wave  packet  is  denoted  by  XI  and  the  band¬ 
width  is  2f. 

It  is  important  to  note  the  distinction  between  6k  and  Ak. 

m.  =  lit  all  -  111 
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The  differential  jk  depends  only  on  the  difference  in  the  magnitudes  of  the 
component  wave  numbers.  The  differential  Ak  depends  on  both  the  difference 
in  the  magnitudes  of  the  component  wave  numbers  and  their  directions.  The 

♦  4 

directions  of  both  k  and  6k  with  respect  to  the  positive  x-axis  is  Y.  The 

«* 

direction  of  Ak  with  respect  to  the  positive  x-axls  is  denoted  by  8.  It  can 
be  shown  (Breeding,  1978)  that 

<S4L  =  A4L  cot.  <|> 

where 

*  =  e- y 

The  velocities  of  the  wavelets  and  packets  relative  to  the  current  are 
determined  by  holding  the  phases  constant,  respectively,  in  (2).  The  wave¬ 
let  (phase)  velocity  is  defined  by 


V-  =  4  I 


4l.  '■*. 


where  e^  is  a  unit  vector  in  the  direction  of  k.  The  velocity  of  the  wave 
packets  is  called  the  geometric  group  velocity  (Breeding,  1978),  and  it  is 
given  by 


L.  i  9(0  _  .  i  wi 

G  *  e- 


Vws. 


where  e  is  a  unit  vector  in  the  direction  of  Ak. 
m 


When  (1)  is  solved  for  r'  and  the  result  is  substituted  into  (2)  the 
displacement  is  expressed  by  A  +  f 


SUlttLCWCUI.  XO  L CO0CU  UY  _ _  .  dk  ^  1 
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From  (10)  it  is  seen  that  the  frequencies  in  the  fixed  coordinate  system  are 
defined  by 

W  s  0)1  +  X- u- 

t  =  qr'  + 


Despite  the  difference  in  frequencies,  the  wavelengths  ^ ■  27T/k  and  L  ■  27T/m 


are  the  same  in  both  the  xy  and  x'y '  coordinate  systems.  The  velocities  of 
the  wavelets  and  wave  packets  in  the  xy  coordinate  system  are  found  by  taking 
the  time  derivative  of  (1)  and  substituting,  respectively,  (8)  and  (9)  to 
obtain 

\r  s  /\r*  -v  u. 

G  =  G  +  u. 

where  v'-<»>Vk  and  G1  ■  (OCjJ'Ok)  cos  $.  If  the  current  vanishes  (13)  and  (14) 
with  definitions  (8)  and  (9)  are  the  same  as  obtained  by  Breeding  (1978)  for 
the  velocities  of  the  wavelets  and  wave  packets,  respectively. 

3.  Snell's  laws  for  hydrons  for  refraction  due  to  bottom  topography 
and  currents 

Wave  refraction  is  considered  for  the  idealized  case  where  a  current  is 
directed  parallel  to  a  parallel  set  of  water  depth  contours.  Both  the  water 
depth  and  the  magnitude  of  the  current  vary  in  the  lateral  direction.  The 
refraction  of  the  wavelets  is  illustrated  in  Figure  1.  A  discontinuity  in 
both  the  water  depth  and  current  magnitude  are  indicated.  The  subscript  i 
denotes  the  incident  wave  while  the  subscript  t  depicts  the  transmitted  (re¬ 
fracted)  wave. 

The  wavelet  crest  must  be  continuous  across  the  boundary.  Therefore, 
as  the  Incident  wave  crest  advances  one  wavelength  the  refracted  wave  crest 
also  advances  one  wavelength.  With  respect  to  the  fixed  coordinate  system 
the  time  taken  for  a  crest  to  advance  one  wavelength  is  the  same  as  the 
incident  wave  period  on  both  sides  of  the  boundary.  However,  relative  to  the 
currents  the  wave  period  changes  as  the  waves  cross  the  boundary. 

The  angle  of  refraction  is  determined  from  the  condition  that  the  com- 


ponent  of  Che  wave  number  parallel  Co  Che  boundary  is  Che  same  on  boch  sides 
of  Che  boundary.  The  resulc  is  Snell's  law 

JjL  A wv\.N  =  omA  (15) 

The  angle  Y  is  Che  angle  Che  cresC  makes  wich  Che  boundary  or  alCemaCively 
Che  angle  between  Che  orChogonal  Co  Che  cresC  and  Che  normal  Co  Che  boundary. 

Since  remains  consCanC  in  Che  fixed  coordinace  system  (15)  can  be  divided 


by  <0  to  obCain 


where  wich  appropriace  subscripts 


V  =  V  + 

The  phase  speed  v'  is  defined  relative  Co  Che  current,  and  u^  is  the  com- 
ponent  of  the  current  in  Che  direction  of  k.  The  component  u^  is  given  by 

U.«  =  u  WY 


(lb) 


ui) 


Cl*) 


The  refraction  law  can  be  stated 


+  ULtAi^Yt  *  Vj  +  mLAi*Of‘L 


(W 


If  u^  ■  0  this  resulc  reduces  to  the  one  obtained  by  Johnson  (1947). 

The  refraction  law  for  hydrons  is  obtained  in  the  same  manner  as  is  the 
wavelet  refraction  law.  A  crest  determined  by  Che  interference  maximum  of 
Che  hydron  (alternatively,  the  nodal  line  between  wave  packets  can  be  con¬ 
sidered)  must  be  continuous  across  a  boundary.  Snell's  law  for  the  hydron  is 


yvA.  avw.  8  = 

where  6  is  Che  direction  of  m  wich  respect  to  Che  normal  Co  Che  boundary. 

In  Che  fixed  coordinate  system  7  is  constant  on  both  sides  of  the  boundary. 
When  (20)  is  divided  through  by  U  Che  resulc  is 

—  ■  -  <ho 

6* 


(3.0) 


where  the  differencial  forms  of  (4)  and  (6)  have  been  used.  Further 


6*-  <3  +  u, 


(**-) 


where  G'  is  the  speed  of  the  hydron  relative  to  the  current  and  u  ■  u  sin  e 

IQ 

is  the  component  of  the  current  in  the  direction  of  m.  The  hydron  refraction 
lav,  which  is  the  counterpart  of  (19),  can  be  stated 


Gt'  +  ”  G-J  +  8i 


(13) 


where  G'  is  defined  by  (9) . 

In  the  absence  of  any  currents  (23)  reduces  to  the  results  obtained  by 
Stoneley  (1935)  and  Breeding  (1978).  Stoneley  used  the  result  to  determine 
the  direction  of  amplitude  (group)  fronts  along  monochromatic  trajectories, 
whereas  Breeding  used  the  result  to  define  wave  packet  trajectories. 


4.  Ray  theory  and  Hamilton's  equations 


There  is  a  very  useful  analogy  between  ray  theory  and  particle  mechanics. 
Synge  (1962)  has  used  the  dispersion  relation  for  water  waves  to  derive  an 
equation  which  is  an  analogue  of  the  Hamllton-Jacobi  equation  used  to  solve 
mechanical  problems.  Goldstein  (1950),  Landau  and  Lifshitz  (1959),  and 
Lindsay  (1960)  discuss  the  similarity  between  the  eikonal  equation  for  rays 
and  the  Hamilton-Jacobi  equation.  In  what  follows  the  ray  theory  treatment 
presented  in  articles  66  and  67  of  Landau  and  Lifshitz  (1959)  for  nondisper- 
slve  waves  will  be  extended  to  dispersive  waves. 


A  wave  packet  has  two  phases.  The  phased  of  the  wavelets  can  be 


stated 


If  =  —  cot 

The  wave  number  k  and  frequency  fa)  can  be  defined  in  terms  of  y. 

i  =  vy 

=  “II 


CW) 


(8-5} 

(**>) 


By  definition  k  ■  It  +  k  *  Thus,  substitutions  for  k  ,  k  ,  and  CJ 


yield 


(§f*  tyt-±\<kT-  o 


Additionally,  the  wave  packet  has  the  phase  T  where 

T  =•  avvxX  +  -  O’  t 

Further 

A-=  VT 


t>t 


0 

c 


Since  m 


2  2  2  2 

m  +  m  *  (O’  /G  )  it  is  seen  that 
x  y 


*A8|f-^(£r=o 


a 


Equations  (27)  and  (31)  are  first  order  partial  differential  equations 
which  are  similar  to  the  Hamilton-Jacobi  equation.  In  mechanics  S  is  the 
action  of  a  particle.  The  momentum  p  -  VS.  Hamilton's  function  H  is  the 
energy  of  a  particle  and  H  -  -3S/3t.  The  solution  of  the  Hamilton-Jacobi 
equation  is  equivalent  to  solving  Hamilton's  equations 

#=  -JH  = 

^  =  V*-H  *  ■§£ 

where  in  (33)  V  is  the  particle  velocity.  In  ray  theory  the  phases  ^  and  X 
are  analogous  to  S.  The  wave  numbers  ic  and  m  are  analogous  to  and  W  and  U* 
are  analogous  to  H.  Hamilton's  ray  equations  can  be  expressed  as 

-  a 

In  (34)  k  and  r  are  independent  variables  whereas  in  (35)  m  and  r  are 


n.  r 

* 


Independent  variables.  In  the  absence  of  currents  W  is  constant  along  a  ray 
(Landau  and  Llfshitz,  1959).  Equation  (36)  is  the  velocity  of  the  moving 


interference  pattern  which  constitutes  the  wave  packet.  If  there  are 
currents  (36)  becomes  the  same  as  (9) . 

Equation  (34)  can  be  derived  on  the  assumption  that  there  is  conserva¬ 
tion  of  the  wave  crests  (Phillips,  1977).  The  partial  derivative  of  (25)  is 
taken  with  respect  to  time  and  (26)  is  used.  Since  k  and  r  are  independent 
variables  dk/dt  ■  3k/ 9t.  The  result  is  (34).  The  one-dimensional  form  of 
this  equation  was  derived  by  Rossby  (1945).  In  a  similar  fashion  (29)  and 
(30)  can  be  combined  to  obtain  (35) . 

5.  Ray  curvature  formulas  for  wavelets,  packets,  and  rays 

Ray  curvature  expressions  will  be  developed  based  on  the  method  pre¬ 
sented  in  articles  66  and  67  by  Landau  and  Lifshitz  (1959).  It  is  assumed 
that  v,  G,  and  u  are  functions  of  x,y.  In  addition  u  «  v,  u  <c  G,  and 
terms  higher  than  first  order  in  u  are  neglected. 


5.1  Ray  curvature  for  wavelets.  Equation  (11)  is  substituted  into  (34) 


to  obtain 


mL  ,  -  i-W  -  ^  (>•£) 


131) 


where  v'  *  |*»' / k.  It  is  convenient  to  employ  the  vector  identity 

=  t.V*  +  4LXtVKa)  +  +  U-X(.vx£)  (3£) 


The  last  two  terms  on  the  right  hand  side  of  (38)  are  zero  since  k  and  r  are 
independent  variables.  Therefore 

-£•  VvL-i-X  (VXUO  I 

An  alternative  expression  can  be  derived  for  dk/dt.  Let  the  unit 


W) 


vector  e,  ■  k/k.  Then  -»  a 

%  WL  .  .  iefc  .  %  Ml 

V5T 


WO) 


s 


when  (39)  and  (40)  are  combined  and  Che  result  simplified  it  is  found  that 


ft  ^  »  •»  *  I 

=  -  VV  -  V7tL  "  ~  IF 


JLt 

ft  A 

Since  e^e^  ■  1  it  follows  that 


K-  =  o 


Thus  e^  and  de^/dt  are  perpendicular.  The  unit  vector  n^  is  taken  tangent  to 


the  wavelet  crest,  i.e.,  perpendicular  to  e.  .  Then 

A  * 


15-  =  =  -  A-4L*  [Vv‘  +  V  V&  +  Cvxu.)] 


where  Y  denotes  the  direction  of  k  with  respect  to  the  positive  x-axis. 
Equation  (43)  simplifies  to 


a*  V  9V*  av'  j.  a* 

S£  -  ax  -  a*  +  v  IT  ■ 


I 


If  v'  is  constant  (44)  reduces  to  Zermelo's  result  reported  by  Arthur  (1950). 

The  component  of  the  current  in  the  direction  of  k,  denoted  by  u^,  is 
given  by 


A  -* 


Uji  =  e„  •  u. 


4. 

It  can  be  shown  that  (44)  is  equivalent  to 


( 


■fj  -  cv+  f^Cv’  +  u*)  c 

This  result  was  obtained  by  Arthur  (1950). 

The  arc  length  along  the  trajectory  of  an  orthogonal  to  a  wavelet  crest 
is  defined  by  ds^  •  (v*  +  u^)  dt  S  v*  dt.  As  a  result  (46)  becomes 

ktv'+^  -  < 

This  expression  defines  the  ray  curvature  of  the  wavelets,  and  is  a  differ¬ 
ential  form  of  Snell's  law  (19).  If  there  are  no  currents  (47)  reduces  to 
the  expression  derived  by  Hunk  and  Arthur  (1952)  and  Arthur  et  al.  (1952). 


hydrons  is  derived  in  the  same  manner  as  is  the  ray  curvature  expression  for 
the  wavelets.  It  is  found  that 


M  =  0  (G*  +  U^vO  -  I^CG'  +  u^ 

where  6,  the  direction  of  the  hydron,  is  in  the  direction  of  m  with  respect 

to  the  positive  x-axis,  and  u  is  the  component  of  the  current  in  the  di- 

m 

rection  of  m.  The  ray  curvature  of  the  hydron  is  defined  by 

^  =  G1  lx  +  +  CtS) 

where  ds  is  an  arc  length  along  the  trajectory  of  an  orthogonal  to  a  hydron 
in 

crest  (interference  maximum  of  wave  packet) .  This  result  is  comparable  to 
Snell's  law  (23).  In  the  absence  of  currents  (49)  reduces  to  the  result 
presented  by  Breeding  (1981). 


5.3  Ray  curvature  for  rays.  Equation  (14)  can  be  written 

0  =  G'  +  S. 

Consider  the  derivative  d(mG)/dt  where  (50)  is  used. 

=  g1  Asjp  +  ^  +  5*  ^ 

For  a  wave  packet,  in  place  of  (39)  it  follows  that 

=  -  /vv^VG1  -/WwXCVXtx) 

For  steady  state  3G'/3t  ■  0  and  du/dt  *  0.  Therefore,  along  a  ray 

4  =  !•«. 

<Lt 


(.50) 

C511 

(53i) 

C55) 

C5*tt 


When  (50)  is  substituted  into  (54) ,  and  only  the  first  order  term  in  u  is 
retained,  it  is  found  that 


£SL  =  vi  (55) 

<$Lt  /vw 

The  substitution  of  (52),  (53),  and  (55)  into  (51)  yields 

—  C/w\»G)  =  -  g'/vwVG1  —  g'^x  (VX&)  +  (G  «VG  )  +  u.  C5fa) 

J.t 

Since  u  «  G,  and  by  neglecting  terms  higher  than  first  order  in  u,  an 
approximation  to  (56)  is  obtained 

^  C/wv.G')  -  -  Q)  aa.76*  -  (.7  XU.)  +  yw^^O*  VG  )  +  VA.  CS7) 


where  the  unit  vector  e 


G/G. 


Also 


^_.^/vwG)  =.  ^_.(/va.6)  + 


/w\G 


A.t 


05S) 


where  e^  and  de^/dt  are  perpendicular.  A  unit  vector  normal  to  e^  will  be 
denoted  by  n  .  Equations  (57)  and  (58)  can  be  combined  to  obtain 

TET  =  -  VG’  -  ^  +  CM) 

The  third  term  on  the  right  hand  side  of  the  equation  can  be  neglected  since 
u  «  G  and  since  the  wave  properties  are  assumed  to  change  little  over  dis¬ 
tances  of  one  wavelength. 

The  direction  of  the  ray,  l.e.,  the  direction  of  e^  with  respect  to  the 
positive  x-axis  is  denoted  by  P .  Then 

=  /*n.*T§T  =  “  An.*  VG*  -  /*»„,•  X  (Vx  U.)]  CfeO) 

♦ 

For  surface  water  waves  only  the  vertical  component  X,  of  the  vortlclty  affects 
the  ray  trajectory.  The  z-component  of  (7  X  u)  is  given  by 

?•*($*-«»)  (U) 


Simplification  of  (60)  leads  to 

.1 


M.  -  ju^p—  -  c«.p2i  +  42*  _  jas 
jet  ~  3x  1  ay  r  ax  a* 


CU) 


Further,  since  ds^s;  G'  dt,  where  ds^  is  an  element  of  arc  length  along  the 
ray,  it  is  found  that 


£P  _  I 


p 


SSl 

3X 


( 


This  is  the  ray  curvature  expression  for  the  rays.  If  the  currents  vanish 
(63)  reduces  to  the  result  obtained  by  Breeding  (1981). 


6.  Discussion 

Equation  (63)  defines  the  trajectories  of  rays,  and  yields  the  paths  of 
constructive  interference  of  waves  (wave  packets)  accounting  for  the  refrac¬ 
tion  effect  of  bottom  topography  and  currents.  At  each  point  along  a  ray 
the  packet  direction  0  is  determined  by  (49)  and  the  wavelet  direction  Y  is 
determined  by  (47) .  The  behavior  of  the  waves  is  more  complicated  if  there 
is  dispersion  than  if  there  is  not. 

The  phase  speed  v  and  conventional  group  speed  U  -  3GJ/0k  of  a  surface 
gravity  water  wave  for  water  of  arbitrary  depth  h  are  given  by  (Lamb,  1932) 


«>H 

u- tO  ♦*£&*) 

(V 

If  there  are  no  currents  the  velocities  U  and  v  are  in  the  same  direction 
(Phillips,  1977;  Breeding,  1978).  If  the  waves  are  carried  by  a  current  then 
U’  and  v'  are  in  the  same  direction  but  U  and  v  are  not. 

Under  certain  conditions  (64)  and  (65)  can  be  simplified.  If 


tanh  kh  »  1,  a  condition  which  defines  deep  water,  then 


Shallow  water  can  be  defined  by  the  condition  tanh  kh  -  kh.  In  this  event 

\J  -  V  -  «>S) 

In  shallow  water  the  waves  are  not  dispersive.  For  all  other  water  depths 
the  waves  are  dispersive. 

6.1  No  dispersion .  Without  dispersion  (14)  becomes 

G  =  +  u  W) 

The  ray  curvature  expressions  (47)  and  (49)  are  equivalent,  i.e.,  ©•)(  . 

However,  the  rays  are  not  normal  to  the  wave  fronts,  i.e.,  P *  Y .  Arthur 
(1950)  considered  the  refraction  of  shallow  water  waves  moving  in  a  current 
over  a  topography  with  variable  water  depth.  In  the  event  there  are  no 
currents  the  ray  curvature  expressions  (47),  (49),  and  (63)  are  identical. 

Then  P  ■  0  -  Y.  If  in  addition  v  is  constant  Y  does  not  change. 

6.2  Dispersion.  When  there  is  dispersion  and  refraction  occurs  due  to  both 
bottom  topography  and  currents  the  directions  P ,  6,  and  Y  are  all  different. 
There  have  been  several  investigations  in  which  refraction  is  considered  due 
to  only  bottom  topography  or  only  currents.  If  there  aren't  any  currents 
the  ray  curvature  expressions  (49)  and  (63)  are  equivalent. 

Johnson  (1947)  used  (19)  to  determine  the  angle  of  refraction  for  deep 
water  waves  arriving  at  an  angle  to  a  linear  current.  In  deep  water  there  is 
no  refraction  due  to  bottom  topography.  Only  the  change  in  the  direction  of 
the  wavelets  was  considered. 

Kenyon  (1971)  considered  the  refraction  of  wave  packets  by  currents  in 
deep  water.  He  used  Johnson's  (1947)  refraction  result  to  obtain  the  di¬ 
rection  of  the  wavelets.  Ray  trajectories  were  determined  using  the  ray 
curvature  equation 


AP  _  s 
■  1? 


no) 


This  result  follows  from  (63)  if  only  the  current  refraction  terms  are  re¬ 
tained  and  G'  ■  U'.  Kenyon's  method  was  also  used  by  Teague  (1974). 

Independently  of  Kenyon  (1971),  Breeding  (1972)  used  what  is  essentially 
the  same  method  to  determine  the  refraction  of  wave  packets  by  bottom  topog¬ 
raphy.  The  wavelet  direction  was  determined  using  Snell's  law  with  phase 
velocity.  The  ray  trajectories  were  determined  using  the  ray  curvature 
equation 
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cm 


where  0  «P.  As  does  (70),  this  equation  follows  from  (63)  if  G  ■  U,  except 
that  the  refraction  terms  involving  bottom  topography  are  retained  instead 
of  the  terms  in  u. 

The  refraction  method  used  by  Kenyon  (1971)  leads  to  a  prediction  of 
different  directions  for  U'  ■  0  -  u  and  v ' ,  whereas  Breeding  (1972)  obtains 
different  directions  for  I?  and  v.  These  results  are  inconsistent  with  the 
definitions  of  the  conventional  group  velocity  and  phase  velocity  for  sur¬ 
face  water  waves.  However,  the  trajectories  obtained  in  this  fashion  can  be 
a  good  approximation  to  those  determined  using  (63) .  This  is  the  case  for 
refraction  due  to  bottom  topography  if  $  ■  0  -H  has  small  or  moderate  val¬ 
ues.  It  is  also  the  case  for  refraction  due  to  currents  if  additionally 
P  St  Q  where  0  is  determined  by  (49) . 

When  refraction  occurs  due  only  to  bottom  topography  (47)  reduces  to 
the  result  obtained  by  Munk  and  Arthur  (1952)  and  Arthur  et  al.  (1952)  for 
waves  moving  with  phase  velocity.  This  ray  curvature  expression  has  been 
widely  used  to  determine  the  trajectories  of  hydrons  on  the  assumption  that 
the  hydrons  follow  the  same  paths  as  monochromatic  waves.  However,  these 
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trajectories  do  not  coincide  with  the  paths  of  constructive  interference 
which  are  obtained  by  the  superposition  of  Individually  refracted  sine  waves 
to  construct  wave  packets. 

Considering  refraction  due  to  bottom  topography.  Breeding  (1978)  used 
(63)  to  determine  wave  packet  trajectories.  The  wavelet  directions  were  de¬ 
termined  using  Snell's  law  with  phase  velocity. 
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